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Abstract. We discuss the extension of a systematic perturbative QCD based analysis to the x — )• 
1 region. After subtracting a number of effects that transcend NLO pQCD evolution, such as 
target mass corrections and large x resummation effects, the remaining power corrections can be 
interpreted as dynamical higher twists. The quantitative outcome of the analysis is dominated by 
the interplay between the value of as in the infrared region and the higher twists. We uncover a 
dual role played by at large Bjorken x that can be used to experimentally extract its value in the 
non-perturbative regime. 
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INTRODUCTION 

"QCD nowadays has a split personality. It embodies 'hard' and 'soft' physics, 
both being hard subjects and the softer the harder." [1] 

QCD's main goal is to describe the structure of hadrons in terms of its fundamental 
degrees of freedom, the quarks and gluons (partons). Hadrons are observable in both the 
initial and final stages of hard processes, while the existence and properties of partons 
are inferred only indirectly. The underlying idea is that because of the smallness of 
the coupling constant at large enough momentum transfer, Q^, or equivalently at short 
distances ^(I/a/Q^), a hard probe sees hadrons as composed of "free" quarks and 
gluons carrying fractions x of the hadron's Light Cone (LC) momentum, with given 
probability distributions. In Deep Inelastic Scattering (DIS) the latter are identified 
with the Parton Distribution Functions (PDFs) [2]. As increasingly shorter distances 
are probed the PDFs shape in x changes due to radiative processes, according to a 
pattern which is calculable with high accuracy within Perturbative QCD (PQCD) [3]. 
Factorization theorems regulate this fundamental property of the theory thereby allowing 
for the short distance behavior to be evaluated using as input the values taken by the 
PDFs at a given initial Q^. The strong coupling regime of QCD remains incalculable; it 
is accessible only with non-perturbative (approximate) methods. 

The separation and yet coexistence of long distance and short distance structure in 
QCD has by now become naturally accepted as part of a "common wisdom framework" 
underlying the interpretation of most experiments, from Deep Inelastic Scattering (DIS) 
to e+e — )■ hadrons, to hadron-hadron scattering. The concept of duality is implicitly 
used to various degrees, meaning, in the most extreme case, that hadronic observables 
are replaced by calculable partonic ones with little more going into the hadronic forma- 
tion phase of the process (from partons to hadrons or vice versa). In a phenomenological 



context, duality purports to study how a number of properties defined from the begin- 
ning of the hard scattering process, are predetermined and persist in the non-perturbative 
stage. A particularly striking realization of duality, known as Bloom Oilman duality [4], 
is observed in DIS, where for large values of Bjorken x > 0.5 (x — Q^/2Mv, M being 
the proton mass and v the energy transfer in the lab system), and for 2^ ~ 5 GeV^, one 
has an invariant mass of < 5 GeV^ (W^ — Q^{l/x — I) +M^), lying mostly in the 
resonance region. While it is impossible to reconstruct the detailed structure of the pro- 
ton's resonances, these remarkably follow the PQCD predictions when averaged over x 
(see e.g. [5] for a review). 

Duality in inclusive scattering is clearly a phenomenological manifestation of the non- 
perturbative to perturbative transition in QCD, whose origins are still largely unknown. 
It is now becoming mandatory to have a fuller understanding of its working, motivated 
on one side by the existence of highly accurate data at large x, and spurred, on the other, 
by the advent of the LHC where previously unexplored regimes in x and will be 
within reach. 

Our studies of the physical origin of duaUty, and of its impact on our understanding 
of the nucleon's structure started a few years ago when we set up a program to quanti- 
tatively extract the scale dependence of the large jc Jlab Hall C data [6]. In the analysis 
performed in Refs.[7, 8] we addressed several effects that have a large impact at large 
X, namely Target Mass Corrections (TMCs), large x resummation effects, and higher 
twists. This work was then completed, and extended to polarized data in Ref.[9]. An 
important point emerged from Refs.[7, 8, 9] that a deeper understanding was needed of 
those aspects unique to the large x perturbative QCD analysis. Our point of view is that 
only after a complete perturbative QCD analysis is performed can one define duality by 
quantitatively establishing whether, and to what extent, this phenomenon is responsible 
for the apparent cancellation of multiparton correlations. This point of view is somewhat 
complementary to the approaches of Refs.[10, 11, 12], which are quark model based and 
focused on symmetry aspects. 

In this contribution we argue that once the range of validity of parton-hadron duality is 
defined quantitatively, a precise PQCD analysis at large jc would open up the possibility 
of extracting the strong coupling constant, as, at low scale. Such an analysis would 
complement the recent extractions using data on the GDH sum rule [13]. It would, 
furthermore, add insight on a recent interpretation put forth in Ref.[14] of the effective 
coupling constant in the strongly interacting/non perturbative regime from light front 
holographic mapping of classical gravity in Anti de-Sittter (AdS) space. 

In order to explain our approach, we first present an overview of the large x data, and 
discuss a few aspects of the evolution mechanism for DIS at large x where two scales 
related to the invariant mass and to the four-momentum transfer, are simultaneously 
present. We reiterate that accurate analyses in this region such as the ones first conducted 
e.g. in [7, 8, 9] are crucial for establishing the interplay of the various components (as, 
multiparton correlations, etc..) of the perturbative to non-perturbative transition regime 
in QCD. We then illustrate the connection between large jc data and as in the infrared 
region, and draw our conclusions. 



ANALYSIS OF LARGE x DATA 



High precision inclusive unpolarized electron-nucleon scattering data on both hydrogen 
and deuterium targets from Jefferson Lab are available to date in the large x, multi-GeV 
regime (see [15] and references therein). Because of the precision of the data one should 
now be able to distinguish among different sources of scaling violations affecting the 
structure functions in addition to standard NLO evolution, 

• Target Mass Corrections (TMC), 

• Large x Resummation Effects (LxR) 

• Nuclear Effects 

• Dynamical Higher Twists (HTs), 

• Impact of Next-to-Next-to-Leading-Order (NNLO) perturbative evolution. 

All of the effects above can be extracted with an associated theoretical error. It is in 
fact well known that their evaluation is model dependent. Recent studies, however, have 
been directed at determining more precisely both the origin and size of the associated 
theoretical error. Recent analyses have been taking into account, so far, some but not all 
of the effects listed above [16]. 



Unpolarized structure function. 

The inclusive DIS cross section of unpolarized electrons off an unpolarized proton is 
written in terms of the two structure functions F2 and Fi, 



dxdy Q^xy 



(1) 



with y — v/£\,ei being the initial electron energy. The structure functions are related by 
the equation 

Fi=F2(l + r2)/(2x(l+i?)), (2) 

where 7^ = AM^x^ /Q^, and R is ratio of the longitudinal to transverse virtual photo- 
absorption cross sections. In QCD, F2 is expanded in series of inverse powers of Q^, 
obtained by ordering the matrix elements in the DIS process by increasing twist T, which 
is equal to their dimension minus spin 

Q^) = ^(x, Q^) + ^ + ^ (^) ^(x, Q^) (l + + ^ (^) 0) 

The first term is the leading twist (LT), with T = 2. The terms of order 1 / Q^~^, T > 4, 
in Eq.(3) are the higher order terms, generally referred to as higher twists [2]. 



Target Mass Corrections 



TMCs are included in Ff^. For > 1 GeV^, TMCs can be taken into account 
through the following expansion [17] 



^1 j^'. 



(4) 



where F2 is the structure function in the absence of TMCs. A more recent analysis [18] 
re-examined TMCs within the coUinear factorization approach of [19] in order to address 
the longstanding question of the unphysical behavior in the threshold region of Eq.(4). 
This originates from the fact that as ;c — )■ 1, one obtains a 2^-dependent threshold, namely 
^2(^ , 2^) for > t„nax = 2/1 + \ + AM^ / Q'^ , therefore rendering F2 undefinable 



as Q varies (see discussion in [20]). In the formalism of Ref.[18], TMCs, appUed to the 
helicity dependent structure functions read 



dm]p{m])F^ 



^1 + 



m' 



(5) 



where Ft = F\. The final quark is assumed to hadronize into a jet of mass mj with a a 
process dependent distribution/smearing function p{inj). In our extraction we take the 
perspective that the evaluation of TMCs is always associated with the evaluation of HTs 
- TMCs should in principle be applied also to HTs - in an inseparable way. Therefore 
we consistently keep terms of ^{l/Q'^) [9, 21], whether in the formalism/prescription 
of Ref.[17] or of Ref.[18]. H{xj Q^), then, represents the "genuine" HT correction that 
involves interactions between the struck parton and the spectators or, formally, multi- 
parton correlation functions. 



Threshold Resummation 



In order to understand the nature of the remaining dependence that cannot be 
described by NLO pQCD evolution, we also include the effect of threshold resummation, 

or Large x Resummation (LxR). LxR effects arise formally from terms containing 
powers of ln(l — z), z being the longitudinal variable in the evolution equations, that 
are present in the Wilson coefficient functions C{z). Below we write schematically how 
the latter relate the parton distributions to e.g. the structure function F2, 

Fi\x,Q^) = ^l,l'dzC{z)q{x/z,Q% (6) 

where we have considered only the non- singlet (NS) contribution to F2 since only 
valence quarks distributions are relevant in our kinematics. The logarithmic terms in 
C{z) become very large at large x, and they need to be resummed to all orders in 
Us. Resummation was first introduced by linking this issue to the definition of the 
correct kinematical variable that determines the phase space for the radiation of gluons 
at large jc. This was found to be = 0^(1 -z)/z, instead of [22, 23]. As a 



result, the argument of the strong coupling constant becomes z-dependent: (Xs{Q^) — )■ 
OCs{Q^i^—z)/z) [24, 25]. In this procedure, however, an ambiguity is introduced, related 
to the need of continuing the value of Us for low values of its argument, i.e. for z — >■ 1 
[26]. Although on one side, the size of this ambiguity could be of the same order of the 
HT corrections and, therefore, a source of theoretical error, on the other by performing 
an accurate analysis such as the one proposed here, one can extract 0C5 for values of the 
scale in the infrared region. We address this point in more detail in the next Section. 

Nuclear Effects 

Theoretical uncertainties in the deuteron are taken routinely into account, and are 
expected to be in sufficient control (see [27] and references therein). Uncertainties arise 
mainly from 

i) Different models of the so called nuclear EMC effect; 

ii) Different deuteron wave functions derived from currently available NN potentials, 
giving rise to different amounts of high momentum components; 

Hi) The interplay between nucleon off-shellness and TMC in nuclei. 

Finally, we did not consider NNLO calculations, these are not expected to alter 

substantially our extraction since, differently from what seen originally in the case of 
F3, these have been proven to give a relatively small contribution to F2. 

Once all of the above effects have been subtracted from the data, and assuming the 
vahdity of the twist expansion, Eq.(3) in this region, one can interpret more reliably any 
remaining discrepancy in terms of HTs. Since we extend ourx-dependent analysis to the 
resonance region we consider the following integrated quantities 

rXram 

r'\{x),Q^)= Fr{x,Q^)dx (7) 

where F2^^ is evaluated using the experimental data in the resonance region. For each ^ 
value: x^^ = / {Q^ + Wl^^-M^), andx„,ax = /{Q^ + W^^-M^), where W^^ and 
W^max delimit the resonance region, and {x) is the average value of x for each kinematics. 
This procedure replaces a strict point by point in x, analysis. 

Typical results from the analysis outlined above are plotted in Figure 1 where we show 
the HT coefficient defined from Eq.(3) as 

R^^ ^ C{x) = [F2ix, Q')/Fi^ix, Q^) - l] 

The error in the figure is from the experimental data. No theoretical uncertainty was 
included. However, our results clearly show that the combined effects of TMCs and LxR 
substantially reduce C{x). We take this as illustrative of the accomplishments one can 
expect from the analysis we suggest in this contribution. Essential features that emerge 
are the interplay between the values of asiM^) and the HTs, the relevance of TMCs, 
and, most importantly, the need to define Us in the infrared region. All of these features 
can affect the central values of the HTs reported in Fig.l. 



H 

in 2.5 



L2.25 
II 

1 2 
c 

J 1.75 
Pi 

1.5 
1.25 
1 

0.75 
0.5 
0.25 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

X 

FIGURE 1. HT coefficients extracted in the resonance region according to the procedure described in 

the text. HT extracted with only the NLO calculation (squares); the effect subtracting TMC (open circles); 
the effect of subtracting both TMC and LxR (triangles). Shown for comparison are the values obtained 
from the coefficient H obtained in Ref.[28] using DIS data and including the effect of TMC. Adapted 
from Ref.[9] 

AT JC ^ 1 

We now discuss in more detail the working of threshold resummation, and its possible 
impact on the analysis of F2 at large x [24]. Starting from NLO, the coefficient, C{z) in 
Eq.(6) is dominated at large x by terms proportional to [<Xs{Q^) ln(l — z)]" which need 
to be resummed in the perturbative series. The physical origin of these terms is in the 
phase space for the contribution of gluons emission to evolution, which become soft as 
X — > 1. A mismatch in the cancellation with the virtual gluons contributions ensues. If, 
however, one carefully evaluates the kinematics for gluon emission at large x within a 
quark-parton model view, one obtains [22] 

q{x,Q^) = ^(^'e^)+/j'^'^/f^..(^'«^)^(^'4), (8) 

where is an arbitrary initial scale, and = 2^(1 — z)/z is the maximum kj_ in the 
virtual photon-quark center of mass system, appearing in the ladder graphs that define 
the leading log result. The resulting phase space is shown for different values in 
Fig. 2. The reduction of the allowed kj^ results in a simultaneous shift in the argument 
of ocs —?■ as{Q^/{l—z)/z), and a cancellation of the as{Q^)ln{l — z) divergence in the 
NLO coefficient function. As a consequence of rescaling the argument of as one has to 
consider its continuation into the infrared region [25, 9]. The left panels of Fig.3 display 
the results for as used in our analysis for different values of Q^. We also show, on 
the right, the extracted value of the effective as from the GDH sum rule. We therefore 
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FIGURE 2. Phase space in the evolution of the NS component. The dotted lines are for kj j^^^ = Q^, 
for g2 = 0.8, 10, 100 GeV2. The full lines represent the upper limit 4,max = for = 10, 100 

GeV^. In this case, one can see a clear reduction of the allowed kj- at large z. 



suggest large x evolution in DIS as yet another way of defining an effective coupling 
constant at low values of the scale. A more quantitative analysis to relate different types 
of measurements, and to study in depth the possible process dependence of as is in 
progress [29]. 

In conclusion, we believe there is a much richer structure to the scale dependence of 
the nucleon's distribution functions that persists behind the apparent cancellation among 
higher twist terms. We started uncovering this structure in the initial work of Refs.[7, 9]. 
While on one side this points at the fact that PQCD provides an essential framework 
for understanding the working of duality, on the other a thorough understanding of the 
lack of final state interactions is still missing. Our analysis opens up the possibility of 
extracting the effective strong coupling constant, as, at low scale from a different process 
thaninRefs.[13, 14]. 

This work is supported by the U.S. Department of Energy grant no. DE-FG02- 
01ER41200. 
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